function dydt=IP_ode_func(t,y,G,P,k)
u=-P.K*y;
d=0.1;
dydt=P.A*y+P.B_u*u+P.B_d*d;
